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Terminal velocities of luminous, early-type SMC stars 

C. J. Evans\ D. J. Lennon\ C. Trundle^'^, S. R. Heap^, D. J. Lindler^ 

ABSTRACT 

Ultraviolet spectra from the Space Telescope Imaging Spectrograph (STIS) 
are used to determine terminal velocities for 11 O and B-type giants and super- 
giants in the Small Magellanic Cloud (SMC) from the Si IV and C IV resonance 
lines. Using archival data from observations with the Goddard High-Resolution 
Spectrograph and the International Ultraviolet Explorer telescope, terminal veloc- 
ities are obtained for a further five B-type supergiants. We discuss the metallicity 
dependence of stellar terminal velocities for supergiants, finding no evidence for 
a significant scaling between Galactic and SMC metallicities for Toff < 30,000 K, 
consistent with the predictions of radiation driven wind theory. A comparison of 
the Voo/v esc ratio between the SMC and Galactic samples, while consistent with 
the above statement, emphasizes that the uncertainties in the distances to galac- 
tic OB-type stars are a serious obstacle to a detailed comparison with theory. 
For the SMC sample there is considerable scatter in Voo/vesc at a given effective 
temperature, perhaps indicative of uncertainties in stellar masses. 

Subject headings: Galaxies: individual: Magellanic Clouds - stars: fundamental 
parameters - stars: winds, outflows - stars: early-type 



1. Introduction 

Observations of individual early-type stars in the Small Magellanic Cloud (SMC) facil- 
itate studies of stellar evolution and radiatively driven winds in a metal-poor environment. 
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Stellar winds are thought to be driven by momentum transfer from metal line absorption; 
thus, one would expect global wind properties to be dependent on metallicity {Z). Partial 
confirmation of this is offered from morphological consideration of SMC ultraviolet spectra 
in which P-Cygni type profiles are weaker (or even absent) when compared with Galactic 
analogs (e.g., Walborn et al., 2000). Previous UV surveys using the International Ultraviolet 
Explorer telescope {lUE; Garmany & Conti, 1985) and the Hubble Space Telescope {HST; 
Walborn et al. 1995) concluded that terminal velocities {voo) for early- type stars are on 
average lower in the SMC than in their Galactic counterparts. This finding has subsequently 
been echoed by Kudritzki & Puis (2000) and Puis, Springmann, & Lennon (2000), albeit 
based on the same observational data. In addition, a theoretical prediction by Leitherer, 
Robert, & Drissen (1992) that Voo scales as is at times referred to in the literature as 
having been confirmed (e.g., Vink, de Koter, & Lamers 2001). 

A detailed analysis of the line driving process for radiatively driven winds, including its 
dependence on metaUicity, was given by Puis et al. (2000). They note that, to first order, 

Voo/Vesc ~ --, T, (1) 

1 — a 

where a is the effective value of the force multiplier parameter, a, defined as the exponent 
of the line strength distribution function. This parameter depends on both metallicity and 
depth (temperature and density) and Puis et al. give explicit examples at various effective 
temperatures and mctallicitics. Provided that we do not enter the regime of weak winds, 
their Figure 27 illustrates that for 0-stars (Tefr = 40, 000 K) the value of a is only weakly 
dependent on metallicity, at least down to one tenth solar. This is confirmed by the calcula- 
tions of Kudritzki (2002), who parameterized the temperature and density variations of the 
force multiplier parameters to construct wind models for 0-type stars to very low metallicity 
(0.0001 Zq). Prom Kudritzki's Table 2 and Figure 9 we see that, at the metaUicity of the 
SMC, one docs not expect any significant change in the value of v^jv^sc for very luminous 
0-type stars (in fact, for the lowest surface gravity models there is even a small increase in 
terminal velocity as metallicity is changed from solar to one-fifth solar). This result does 
not necessarily contradict Leitherer et al. , who note that their scaling is a description of the 
average behavior and that "t'oo niay increase with decreasing Z in certain Teg domains" as a 
result of varying ionization conditions in the wind. For the Kudritzki main-sequence models 
there is indeed a weak trend with Z, which essentially agrees with the Leitherer et al. scaling 
(which for one- fifth solar imphes a change of only 20%, comparable to the uncertainty arising 
from the determination of the surface gravities) . 

In part prompted by the Kudritzki results, here we attempt a more quantitative inves- 
tigation of terminal velocities for supergiants in the SMC; we determine v^o for a total of 16 
targets. Aside from consideration as a whole to study the effects of metallicity, such results 
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are useful individually. Terminal velocities are a necessary "ingredient" for the detailed at- 
mospheric analyses now possible for early-type stars. Theoretical model atmosphere codes 
have undergone significant development over the past decade, most recently with the inclu- 
sion of line-blanketing arising from metallic species e.g., CMFGEN (Hillier & Miller, 1998) 
and FASTWIND (Santolaya-Rey, Puis, & Herrero 1997; Herrero, Puis, & Najarro 2002). One 
of the most salient results from the application of these codes to analyses of 0-type stars 
is the downward revision of stellar temperatures (when compared with previous determi- 
nations that did not include the effects of blanketing), noted for plane-parallel models by 
Hubeny, Heap, & Lanz (1998). This effect has since been more thoroughly explored in both 
the Galaxy and the Magellanic Clouds using extended atmospheres e.g., Martins, Schaerer, 
& Hillier (2002), Crowther et al. (2002), Herrero et al. (2002) and Bouret et al. (2003). 

Detailed line-blanketed methods have currently not been applied widely to B-type stars. 
In combination with high-resolution optical data from the VLT UV- Visual Echelle Spectro- 
graph (UVES), the Voo results presented here for four of our early B-type Space Telescope 
Imaging Spectrograph (STIS) targets are used by Trundle et al. (2004, hereafter TL04) to 
determine intrinsic stellar and wind parameters using FASTWIND. Part of our motivation in 
addressing the dependence of terminal velocities on environment also stems from the TL04 
study. Ultraviolet spectra of four further SMC targets (observed as part of the same STIS 
program) do not display evidence of strong stellar winds and the methods we use here are 
not sufficient for determination of Voo (these additional data will be discussed in a subsequent 
paper). We therefore investigate the scaling relation of with stellar escape velocity (vesc), 
with a view to obtaining estimates of the terminal velocities for these additional targets. We 
note that such scaling relations are also of interest as they are employed in the analysis of 
samples of stars which are too distant or numerous to obtain UV spectra (e.g, Urbaneja et 
al. 2003). 



2. Observations 

2.1. HST-STIS observations 

In the course of our HST General Observer programs (G07437 and G09116, P.I.: 
DJL) we have observed a number of O and B-type SMC stars that display evidence of strong 
stellar winds, for which terminal velocities can be determined. The dwarfs in the young 
cluster NGC 346 (observed as part of G07437) have recently been studied by Bouret et al. 
(2003) and are not considered here. Observational details of our 11 targets are summarized 
in Table 1; all are drawn from the catalog of Azzopardi & Vigneau (AzV; 1975, 1982). 
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Each target was observed using the STIS in echelle mode, with the far-UV MAMA 
(Multi- Anode MicroChannel Array) detector and the E140M grating (centered at 1425 A); 
those from program G09116 were also observed with the near-UV MAMA and the E230M 
grating (centered at 1978 A). All of the observations were made through the 0'.'2 x 0'.'2 
entrance aperture and the effective spectral resolving power of the E140M and E230M grat- 
ings are R = 46,000 and 30,000, respectively. The exposure times for the G09116 data 
were set from consideration of low-resolution lUE spectra. With the E140M grating the 
exposure times range from 54 to 216 minutes (Sk 191 and AzV 18, respectively), with shorter 
exposures for the E230M set-up (ranging from 38 to 75 minutes) . 

The primary CALSTis processing steps are described by Walborn et al. (2000), with the 
difference here that the two-dimensional inter-order background correction performed manu- 
ally in those reductions has now been incorporated into the standard pipeline as SC2dcorr. 
After pipeline reduction the individual echelle orders were extracted and merged to form a 
continuous spectrum. Prior to merger (typically) 25 pixels at both ends of each order are 
clipped; the spectral overlap between orders is generally large enough to accommodate this 
and the final result has a more consistent signal-to-noise ratio than otherwise. The reduced 
echelle spectra cover the wavelength ranges 1150-1700 A (E140M) and 1600-2350 A (E230M). 
For display purposes (and for consistency with Walborn et al. 2000) the spectra are binned 
to 0.25 A per data point. The S/N (on the basis of photon statistics) is typically in the range 
of 20-30 per binned data point. 

In Table 1 we also give H I column densities, A^(H l), toward our targets, found assuming 
the Lya profile arises from a pure damping profile (see e.g., ShuU & Van Steenberg 1985 and 
references therein). To avoid possible contamination by the N V AA1239, 1243 doublet, we 
consider only the blueward wing when deriving the tabulated values. 

2.1.1. Multiplicity of AzV 4? and 216 

Target acquisition images used for the STIS are taken using a 100x100 pixel CCD cam- 
era (covering 5"x5"), with the "longpass" broad-band filter centered at 7230 A (FWHM = 2720 A) . 
Inspection of these images for each of our targets generally show isolated point sources, with 
the exception of a faint companion present in the image of AzV 47. While inspecting the 
acquisition data from the additional targets observed for program 009116,^ it came to our 
attention that the acquisition image for AzV 216 (classified as BlIII by TL04) also displays 



^ Lower luminosity targets were also observed for compilation of spectral libraries and their UV spectra 
are not discussed here. 
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a faint companion. These images are shown in Figure 1. 

The separations are ~3 pixels (0.15") in AzV 47 and ~5 pixels (0.25") in AzV 216 
which, at the distance of the SMC, corresponds to approximately 9000 and 15,000 AU, 
respectively; both companions are relatively faint in comparison to the primary targets. 
Some contamination is possible in the UV spectrum of AzV 47 because the companion will 

not be excluded by the 0"2 x 0"2 aperture. However, Walborn ct al. (2000) noted that 
"the optical and UV spectra of AzV 47 are entirely normal for its type," so it is likely that 
there is not a significant contribution from the companion. A full analysis of the VLT-UVES 
spectrum of AzV 216 is presented by TL04; both Ha and the He I A4471 line show some 
asymmetry, likely originating from the presence of a nearby companion. 

2.2. Additional UV data 

To augment the H ST-STIS spectra, we include four Bl-type stars observed with the 
HST using the Goddard High- Resolution Spectrograph (GHRS) as part of program GO6078 
(P.L: DJL). We also include a high-resolution observation of AzV 78 made with the lUE 
telescope using the large aperture with the short-wavelength prime (SWP) camera (image 
SWP09334, Fitzpatrick & Savage, 1985). The GHRS observations do not cover the Lya 
region, and, although the spectrum of AzV 78 extends far enough into the blue, the signal- 
to-noise in the region of Lya is not sufficiently high to permit measurement of the H I column 
density. These additional targets are also detailed in Table 1. 

3. Determination of terminal velocities: SMC stars 

The stellar winds of B-type stars are much less intense than those observed in 0-type 
stars; terminal velocities (e.g, Prinja, Barlow, & Howarth 1990) and mass-loss rates (e.g., 
Kudritzki et al. 1999) are generally both lower in B-type stars. P Cygni-like profiles are 

seen in the Si IV and C IV resonance lines of three stars (AzV210, 215, and Skl91), and 
weak Si IV emission is seen in AzV 18. Terminal wind velocities were determined for these 
four stars using the Sobolev with Exact Integration (SEl) computer program developed by 
Haser (1995). This program was then employed to find foo-values for the seven 0-type 
STIS targets in Table 1 and for the five archival GBRS/ lUE B-type spectra. The same 
SEI program has been widely used in other recent analyses of early-type stars in the Local 
Group, e.g., Herrero et al. (2001, Milky Way), Bresolin et al. (2002, M31), and Urbaneja et 
al. (2002, M33). 
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The SEI code adopts a standard /3-law to describe the velocity structure and includes a 
turbulent velocity, which increases with radius to a maximum of fturb at v^o (see Haser et al., 
1995 for further details). The C IV AA1548.19, 1550.76 and Si IV AA1393.73, 1402.73 doublets 
are used as diagnostics of v^o- We do not use the N V AA1238.81, 1242.80 doublet in our 
analysis as a consequence of its weakness in early B-type spectra (which constitute the bulk 
of our sample) and because of problems arising from its contamination by Lya absorption. 
The adopted /?-law is essentially a free parameter in the SEI program, to be determined from 
morphological comparisons of the observed and theoretical resonance lines. In the current 
work we are solely interested in the terminal velocity and, though small variations in (5 can 
lead to more aesthetically pleasing line fits, the resulting t'oo-values are still within the quoted 
errors. Therefore, in contrast to other studies (e.g., Bresolin et al., 2002), we do not attempt 
to "fit" the /3-law for our targets and adopt (3—1 (consistent with the results of e.g., Massa 
et al. 2003). 

Our results are given in Table 1, and model fits for four of our B-type STIS targets are 
shown in Figure 2. To correct for the underlying photospheric components of the resonance 
lines in our analysis, we use archival lUE observations of B-type stars with low projected 
rotational velocities (selected from Siebert, 1999) as templates,^ convolved by the appropriate 
V sin i of each target (for estimates see TL04 and Walborn et al. 2000) . The model spectra 
were also convolved with an instrumental broadening function, although given the high 
resolution of these data this effectively has no impact on the derived velocities. 

Terminal velocities for the B-type spectra are determined primarily from consideration 
of the Si IV doublet. However, as shown in the top panel of Figure 2, we are unable to match 
the morphology of the Si IV doublet in AzV 215. The presence of the relatively flat "shelf" 
resembles that seen in k Ori (HD 38771, B0.5 la, Walborn, Nicols-Bohlin, & Panek 1985) 
though in AzV 215 it is more extreme. The C IV proflle for this star is shown in Figure 3 
and is well matched by a, — 1400 km s~^ model. 

In determining Voo for our 0-type spectra we use the C IV proflle as the primary diag- 
nostic; the Si IV doublet is a relatively weak feature in the 0-type stars studied here (see 
Walborn et al. 2000). Correcting for the photospheric component is comphcated in 0-type 
stars by the absence of templates uncontaminated by stellar winds (the same problem is dis- 
cussed by Massa et al. 2003). In Figure 3 we show model fits for two of our 0-type spectra 
(AzV 75 and 95) and also for AzV 215. In the case of AzV 215 the same Voo-value is found 
regardless of whether the photospheric correction is included (using v Ori, HD 36512, Wal- 



^Additional tests with new SMC metallicity tlusty models (R. Ryans, 2003, private communication) 
yield essentially identical terminal velocities. 
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born et al. 1985), i.e., the large terminal velocity of the wind leaves the blue-edge relatively 

free of photospheric contamination. To test if this also holds true for our 0-type stars we 
used the lUE spectrum of the subluminous dwarf BD+75 325 (classified as 05p [Gould, 
Herbig, & Morgan 1957] and used by Herrero et al., 2001, in analysis of their hottest stars) 
as a template for AzV 75 and 95. As demonstrated in Figure 3, the blue edge of the C IV 
profile is again largely independent of the inclusion of a photospheric template, although a 
more visually pleasing fit is achieved to the emission profile. For the remaining 0-type stars 
we simply match the blue edge, without including a photospheric template. 

It is important to note that the determination of the best fit in all our spectra is 
subjective. In the case of saturated lines different combinations of Voo and VtuA can be used 
to match the blue wing of the profiles; in nonsaturated, weaker lines the adopted line-strength 
of each doublet is an additional parameter to be considered. Typical uncertainties in v^^, are 
± 100 km s~^, except in the case of the four B-type stars with v^o < 500km s""*^, where it is 
in the region of ±50 km s~^. 



3.1. Compcirison with previous results 

Terminal velocities for three of our targets were determined by Prinja (1987) from 
comparison of low-dispersion lU E spectra with theoretical profiles. His value for AzV 242 
(1000 km s~^) is within our quoted uncertainty, though his result for AzV 215 (1675 km s~^) 
is slightly larger. More worrying is that his value for Sk 191 (850 km s"^) is twice our result. 
This prompted us to retrieve the relevant lUE spectrum from the archive and it neither 
displays P-Cygni emission nor saturated absorption profiles. From the new higher quality 
STIS data, measurement of the blucward limit of the saturated core of both the C IV and 
Si IV doublets (generally referred to as vuack in the literature) gives a velocity of 450 km s~^, 
proscribing the larger result. 

A more recent study of AzV 242 (Prinja & Crowthcr, 1998) gives Voo = 860km s^^ from 
consideration of observed narrow absorption components (NACs; a method we do not pursue 
here). Albeit at the lower limit, this result is again within our expected uncertainty. 

Garmany & Fitzpatrick (1988) also used low-dispersion lUE spectra to study a sample 
of SMG stars and there is some overlap with our targets (namely AzV 15, 69, 75, and 95). 
They were limited to estimates of the terminal velocities, in effect measuring a v^dge (usually 
defined as where the absorption component of the P-Cygni profile meets the continuum). 
We do not consider such measurements here; the relation of Vedge to Voo is dependent on the 
morphology of the blueward part of the profile, which is primarily determined by the tur- 
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bulent velocity. As one would expect, their velocities are considerably larger (~500 km s ^) 
than those from the new analyses. 

Finally, detailed non-LTE, line-blanketed model atmosphere analyses of AzV 69 and 83 
were presented by Hillier et al. (2003). Prom fitting the absorption and emission line widths 
of selected UV lines, terminal velocities of 1800 km s~^( AzV 69) and 940 km s~^ (AzV 83) 
were obtained from comparisons with model spectra; these results are well within the quoted 
uncertainties when compared with our new SEI results in Table 1. 

4. Determination of terminal velocities: Galactic stars 

To provide a test of our methods we have also used the SEI program to determine termi- 
nal velocities for five Galactic B-type supergiants taken from the Walborn et al. (1985) lUE 
atlas. These stars are drawn from the sample of ~60 early B-type supergiants analysed by 
Hascr (1995), and a comparison of results is given in Table 2. The velocities determined here 
are smaller than those found by Haser, and consequently the adopted turbulent velocities 
are larger (to match the blue edge of each line as described in §3). 

In their analyses of 0-type stars Groenewegen & Lamers (1989), Hascr (1995), and 
Massa et al. (2003) adopted turbulent velocities in the range 0.05-0. 15 foo- Comparable 
results are found for the 0-type stars studied here, with the exception of AzV 83. Hascr 
used similar turbulent velocities for his B-type sample (i.e., 0.05-0. 15 ?;oo)- This highlights a 
potential difference in methods between Haser's and the current study; our adopted turbulent 
velocities are generally larger. Both Bresohn et al. (2002) and Urbaneja et al. (2002) also 
noted that t'turb/'^^oo was larger for their analyses of B-type spectra, compared with previous 
O star results. 

One factor contributing to the larger ratio is simply that the terminal velocities in the 
cooler B-type stars are typically smaller than those in the O stars. Coupled with the fact that 
i'turb does not appear to vary systematically with spectral type this leads to a larger I'turb/'^^oo 
ratio. However, inspection of Table 2 highlights an important issue that arises from our 
Galactic analyses. In all five cases we adopt larger turbulent velocities than Haser to match 
the same observations, using the same program. The situation is most severe in HD 152236, 
where fturb is a factor of 4 larger. In Figure 4 we compare our fit to the HD 152236 Si IV 
doublet with a model calculated using the parameters found by Haser. The larger turbulent 
velocity adopted in our study is clearly more consistent with the blue-edge of each component 
than in the Haser model; note also that the width of the absorption profile is better matched 
by our model. The same features are present when comparing the other four lUE targets 
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with models calculated with Haser's values. The Haser v^o results are generally at the "upper 
limit" of our results, suggesting that his methods were slightly different; i.e., he found the 
maximum velocity consistent with the data, whereas our approach has been to find the range 
of Voo (for which reasonable fits arc found) and then taking the central value. These results 
highlight the dependence of the derived terminal velocity on the adopted turbulence. We 
emphasize that Wturb in the present context is used very much as a fitting parameter, the 
physical interpretation of such large turbulent velocities is extremely uncertain. 

As an aside, we note that the lUE image used for our analysis of HD 91969 is not that 
from the Walborn ct al. atlas (SWP06510). Haser did not specify which image was used in his 
analysis, but his v^o of 1550km s^^ is comparable with previous results from SWP09076, i.e., 
1545 km s~^ (Prinja et al., 1990) and 1500 km s~^ (Lamers, Snow, & Lindholm 1995). The 
wind variabihty of HD 91969 has been discussed recently by Prinja, Massa, & FuUerton (2002) 
and is highlighted by our results; for the SWP09076 spectrum we find Voo — 1500 km s~^, 
cf. 1350 km s~^from SWP06510. 



5. Discussion 

The correlation between terminal and stellar escape velocities was noted by Abbott 
(1978) and subsequent studies (e.g., Groenewegen, Lamers, & Pauldrach, 1989; Howarth & 
Prinja, 1989; Lamers et al. 1995) have revisited this relationship. A rehable, well-constrained 
Voo/vesc relation is an extremely attractive prospect. For instance, if the wind momentum- 
luminosity relationship (WLR, e.g., Kudritzki & Puis 2000) is to be used to determine 
distances to galaxies beyond the Local Group, ultraviolet stellar spectra are generally not 
available (e.g, Urbaneja et al. 2003). In comparison with theoretical spectra, optical data 
can be used to provide an initial estimate of Vesc, thence an estimate of v^o, with some 
iteration to a final solution. Before investigating this relationship for our SMC sample, we 
compare our results with those from Galactic studies. 

5.1. Terminal velocities 

Common belief holds that terminal velocities in the SMC will be lower than in the 
Galaxy because of the decrease in metallicity (e.g., Kudritzki & Puis, 2000). Kudritzki 
& Puis presented a comparison of terminal velocities (as a function of temperature) for 
Galactic and Magellanic Cloud targets. A necessary step in this process was the adoption 
of a suitable calibration to allocate effective temperatures to a given spectral type (in their 
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case that of Humphreys & McElroy, 1984). Over the past 20 years, significant developments 
have been made in model atmosphere techniques, most recently the downward revision of the 
temperature scale discussed earlier owing to the inclusion of line blanketing (e.g., Martins 
et al. 2002). If we are to compare our foo-values with Galactic results in this manner, it is 
clear that any differences in temperature scales should be minimized. 

Temperatures for the SMC targets are taken from line-blankctcd, non-LTE analyses by 
Hillier et al. (2003; AzV 69, 83), TL04 (Sk 191, AzV 18, 210, 215), and Trundle et al. (2004, 
in preparation; AzV 78, 242, 264); note that for consistency with the analysis of AzV 69 
and 83 by Hillier et al. wc adopt their Woo-valucs. Temperatures for the remaining targets 
(AzV 15, 75, 95, 327)'^ are estimated from interpolation between the line-blanketed results 
of Crowther et al. (2002) and Martins et al. (2002). Our results are supplemented by those 
from similar analyses of SMC targets by Crowther et al. (2002; AzV 232) and Evans et al. 
(2004; AzV 70, 235, 372, 456, 469, 488). Terminal velocities for this expanded sample are 
plotted as a function of stellar effective temperature (Teff) in Figures 5 and 6. 

In Figure 5 we compare our results with the data of Lamers et al. (1995). Their 0-type 
temperatures were taken from the unpublished Garmany catalog and Chlebowski & Garmany 
(1991), in which non-LTE results were used to define a temperature scale. Temperatures for 
their B-typc targets were taken from the LTE calibrations of Schmidt-Kaler (1982). From 
qualitative inspection of Figure 5 it appears as though the SMC values are generally lower 
than the Galactic values, as previously assumed in the literature. 

In Figure 6 we show our results together with recent results from detailed analysis 
of individual Galactic targets, i.e., independent of any external temperature calibrations. 
Results for 0-type supergiants (and giants) are taken from the non-LTE, line-blanketed 
studies of Herrero et al. (2002) and Repolust, Puis, & Herrero (2004). TL04 presents such 
analyses for B-type SMC supergiants, but there arc no contemporary results for Galactic 
B stars. In Figure 6 we include results from McErlean et al. (1999), which employed non- 
LTE methods to analyze Galactic B-type stars, but which predate the availability of line- 
blanketed atmosphere codes; values are taken from Haser (1995) or Howarth et al. (1997). 
As discussed in TL04, for the early B-type spectra the inclusion of hne- blanketing has a 
relatively minor effect on the temperature scale (cf. the results of Dufton et al. 2000), i.e., 
McErlean's temperatures are preferable to those assumed from Schmidt-Kaler. 

In contrast to the comparison with the Lamers et al. sample, there is no obvious 
qualitative difference between the Galactic and SMC results for Teff < 30, 000 K. Linear 



■'We omit AzV 80 from this discussion since the exact nature of the peculiar "nfp" spectral class is unclear 
(Walborn et al. 2000). 
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fits to these results (excluding the Galactic A-type supergiants) give slightly different re- 
lations for the Galactic and SMC samples. However, because of the scatter of the results, 
these differences are not statistically significant and therefore we do not undertake a rigorous 
quantitative analysis. 

At hotter temperatures the situation is somewhat more ambiguous; our motivations 

for the current study were primarily concerned with B-type supergiants and the sampling 
of targets with T^s > 30, 000 K is relatively sparse. In Figure 6 we distinguish between 
luminosity class 1 and III objects. With the exception of the low velocity for HD 191423, 
the Galactic giants appear to have generally larger tioo-values than the Galactic supergiants. 
A more objective comment is that Voo for the SMC giants is lower than for the Galactic 
targets (again excepting HD 191423), although from such a limited sample such differences 
clearly cannot be attributed to the change in metallicity. We are similarly limited by our 
sample of more luminous 0-type stars. Both AzV83 and 232 have low terminal velocities 
for their spectral types (noted by Walborn et al. 2000) and appear distinctly different from 
the Galactic supergiants. However AzV83 and 232 are both "laf" type spectra with dense 
winds (e.g., Hillier et al. 2003); foo- values for similarly extreme stars in the LMC are also 
relatively low, e.g., Sk — 66°169 (Crowther et al. 2002, Massa et al. 2003). The Galactic 
supergiants in Figure 6 are less extreme objects and therefore do not permit a sensible direct 
comparison of 0-type supergiant terminal velocities at the current time. 



5.2. Stellar escape velocities 

We now consider the escape velocities for SMC stars with published intrinsic stellar 
parameters (see Table 3). As is usual in such studies, we define the effective escape velocity 
(vesc) as the gravitational escape velocity, corrected for the contribution to the radiative 
acceleration caused by Thomson scattering by free electrons. Therefore, 



l2 GM^(l-r) 

where F is the ratio of the radiative acceleration from electron scattering to the stellar gravity 
given by 

AirGM.c ^ ' 

in which (Tg is the electron scattering opacity, given by n^o^jp^ where n^, is the electron 
number density, ctt is the Thomson cross-section, and p is the density of the gas. From 
the TL04 fastwind model for AzV 215, = 0.31 (at r = 2/3). This value is consistent 
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with previous values in early- type stars (e.g., using equation (3) and their values for F, — 
0.28-0.32 in the Galactic study of Lamers et al., 1995) and is used here in our calculations. 

Stellar radii, masses, and luminosities are taken from the relevant papers; note, however, 
that the SMC distance modulus adopted by Hillier et al. was 19.1 (cf. 18.9 in the other 
studies) which will affect the luminosity of the final atmospheric models, hence the radii and 
masses and ultimately the escape velocity (of order 10%). 

In Figure 7 we compare v^jv^^ for the SMC stars in Table 3 with the results from 
Lamers et al. (1995). As discussed before in §5.1, the temperature scales are not directly 
comparable but both the scatter and the mean value of Voo/f^esc is roughly consistent with 
that of the Lamers sample. From a compilation of published results, Kudritzki & Puis (2000) 
found Woo = 2.65 Vesc (Tcs > 21,000 K) quoting an accuracy of ~20%. The mean from our 
sample (excluding AzV 18 and 210 as they are cooler than 21,000 K) is 2.36; AzV83 is again 
the most distant outlier. 

In Figure 8 our results are compared with those for Galactic 0-type giants/supergiants 
(Hcrrero et al. 2002; Repolust et al. 2004) and B and A-type supergiants (Kudritzki et al. 
1999). The results for HD 207198 and 209975 (Repolust ct al. 2004; v^jv^^^ = 6.09 and 
5.62, respectively) are excluded from the figure. These two stars are listed by Markova et 
al. (2004) as being members of the Cep 0B2 association, for which they adopt a distance of 
0.85 kpc. However, Hipparcos results by de Zeeuw et al. (1999) support a smaller distance of 
0.62 kpc. The region is also extended over several degrees, so there will hkely be a significant 
depth effect within the association, leading to uncertain relative distances. From equation 
(2) overestimates of their luminosities will lead to underestimates of Vesc- Smaller distances 
for these two stars would bring them into better agreement with the other results in Figure 
8. The mean Voo/vesc of the remaining Galactic targets (Teg > 21,000 K) is 3.07, more in 
keeping with Abbot's (1978) suggestion that v^o ~ 3fesc than the Kudritzki & Puis result. 

Uncertainties in the distances to the Galactic targets may contribute to these differing 
results and also to the greater spread of foo/'^csc (cf. Lamers et al.). However, it is worth 
noting that the targets in the Herrero et al. (2002) study are all (thought to be) members of 
Cyg 0B2; aside from any systematic offset, the spread of results for their seven targets should 
be minimized and yet a wide range of Voo/'^esc is still found. Unfortunately, from current 
methods there will always be a large scatter associated with Voo/ Vesc, primarily a result of 
the uncertainty in logs' of (at least) ±0.1 dex. This corresponds to a large uncertainty in 
the spectroscopic mass, which in the case of, e.g., AzV 215, yields results for Vesc of ±20%. 

In addition to its use for targets without ultraviolet information, understanding the 
scaling of Voo with Vesc is useful for targets such as AzV 22 and 362 (see TL04 for analysis). 
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The optical spectra of these two stars display evidence of a stellar wind, but their ultraviolet 

spectra (also observed as part of G09116) are devoid of strong wind signatures. Similarly, 
infilling of the Ha profiles is observed in AzV 104 and 216 indicative of some, albeit weak, 
stellar wind. Since we are limited in the current work by a relatively small sample, TL04 
used the Kudritzki & Puis formula to obtain v^o estimates for inclusion in fastwind. 

In general, the results in Figure 8 display no evidence for a significantly different t'oo/^csc 
ratio in the SMC in comparison with that for Galactic stars. We therefore conclude that in 
the case of supergiants with strong winds, our results are in good agreement with predictions 
of theory. 

Figures 7 and 8 also prompt further discussion in the context of the so-called "bistability 
jump" near T^s = 21,000 K (Lamers et al. 1995). From Figure 7, there appears to be a clear 
step between the high (hotter) and low (cooler) ionization regimes, but such a well-defined 
edge is largely an artifact of the temperature assignment for the spectral bins. The Lamers 
et al. sample includes seven Bl-type stars for which the same temperature (20,800 K) is 
used; in reality, the temperature for a given spectral type encompasses a range of values, 
e.g., Skl91 and AzV 210 are both B1.5 la-type spectra, yet their temperatures differ by 
2000 K (TL04). Further comphcations can arise from the effects of hne-blanketing, i.e., 
the temperature derived for a star with a strong wind can be lower than that found for a 
later spectral type with a much weaker wind. In these cases the temperature scale is not 
necessarily monotonic (e.g., TL04, Evans et al. 2004) and the temperature range for a given 
spectral type can overlap with that of another. In Figure 8 there is still evidence for a lower 
foo/fesc ratio at cooler temperatures, but, as a consequence of analysing individual stars 
rather than employing an effective temperature-spectral type calibration, the boundary is 
somewhat "blurred." Additionally, we note that the median Voo/i'esc for those SMC stars 
significantly hotter than this region (Teff > 24,000 K) is 2.63, consistent with the Kudritzki 
& Puis relation. 



6. Conclusions 

From high-quality ultraviolet spectra we have used the SET method to determine termi- 
nal velocities for nine B-type supergiants and seven 0-type stars in the SMC. One important 
conclusion from the current work is that, if one is solely interested in terminal velocities, 
many of the parameters in the SEI code are essentially cosmetic fitting parameters to which 
other authors have attributed undue significance in previous studies. 

We compare our results with those from recent Galactic studies, adopting temperatures 
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from tailored atmospheric analyses, independent of external calibrations. To our knowledge, 
this is the first time such a thorough comparison has been made, giving a clearer view of 
the scaling of terminal velocities. We find no evidence for lower terminal velocities in the 
SMC than in the Galaxy for T^s < 30,000 K (corresponding to spectral types later than 09). 
Similarly, in agreement with the predictions of theory for luminous supergiants, we find no 
evidence for a significantly different scaling of v^o with Vesc in the SMC in comparison with 
that in the Galaxy. We stress that the situation for giants and dwarfs remains untested. 

From our limited sample it appears that Voo for 0-type giants may be lower in the SMC 
than in the Galaxy, although both this and the situation for 0-typc supergiants warrants 
further investigation. Further results from contemporary model atmosphere analyses of O- 
type stars in the Galaxy and the Magellanic Clouds would greatly help attempts to delineate 
the scaling of v^o with both temperature and Vgsc- 

Finally, given the uncertain distances, it would appear that calibrating the WLR using 
Galactic OB-type stars is not an attractive prospect at the current time. 
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Table 1. Observational parameters of SMC targets 
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References. — (1) Azzopardi & Vigneau 1975, 1982; (2) Garmany et al. 1987; (3) Massey 2002; (4) Walborn et al. 2000 

Note. — Stellar identifications arc those of Sandulcak (1968, Sk) and Azzopardi & Vigneau (1975, 1982, AzV). Classifications for the O-type spectra are 
from Walborn et al. (2000); B-type classifications are from Lennon (1997), excepting AzV 483 (Lennon, 1999) and AzV 78 (from re-inspection of the blue 
optical data the spectrum is now classified as Bl Ia-|-). Radial velocities are found from the blue optical data from Lennon (1997) and an unpublished echelle 
spectrum of AzV 483 (using CASPEC at the ESQ 3.6-m telescope). For AzV 18, 78, 483 and Sk 191 (i.e., Voo < 500) the uncertainty in Voo is ±50 km s'^and 
for the remainder of the sample this increases to ±100 km s~^. 
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Table 2. Terminal velocities determined for selected Galactic B-type supergiants 
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Note. — All tabulated velocities are given in km s , typical uncertainties of v^o are 
±100 km s-^ 
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Table 3. Escape velocities for O and early B-type SMC stars 



Star 
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Ref. 
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Note. — Stellar temperatures have been determined 
using non-LTE, line-blanketed model atmospheres. The 
values for AzV 83 and 69 are flagged as uncertain as 
a result of the larger SMC distance adopted in their 
analysis compared with the other studies. 

References. — (1) Crowther et al. 2002; (2) Evans et 
al. 2004; (3) HiUier et al. 2003; (4) Trundle et al. 2004, 
(5) Trundle et al. 2004, in preparation. 
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Fig. 1. — STIS acquisition images for AzV 47 (left) and AzV 216 (right). The regions shown 
correspond to approximately 3'.'8x2'.'9. 
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Fig. 2. — Model fits to the Si IV doublet for STIS B-type spectra. We are unable to match 
the morphology of AzV 215, though the same Voo successfully matches the C IV profile (see 
Figure 3). 
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Fig. 3.— Model fits to the C IV doublet for AzV 75 [05 III(f+)], AzV 95 [07 III((f))] 
and AzV 215 [BNO la]. Models were calculated without consideration of the photospheric 
contribution (left panel) and with inclusion of a photospheric template (right panel). The 
derived terminal velocities are identical. 
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Fig. 4. — Si IV doublet of HD152236 {solid lines) compared with model profiles from the 

SEI code {dot-dashed line). The top panel shows a model calculated for Vco = 500 km s~^, 
Vturh = 50 km s~^and /? = 1.5 (Haser, 1995). The lower panel shows the model from the 
current study: v^o = 450 km s^^, fturb = 175 km s^^'^and P = 1.0. All profiles are shown in 
velocity space, relative to the blue component of the doublet, i.e., 1393.7 A. 
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Fig. 5. — Stellar terminal velocities, function of effective temperature for the 

current SMC sample {solid symbols) and for the Galactic targets of Lamers et al. (1995, 
open symbols). Squares are those stars for which the temperatures have been interpolated 
from published results. 
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Fig. 6. — Terminal velococity function of Teg for the current SMC sample {solid 

sym,hols) and for published analyses of Galactic targets {open symbols, see text for sources). 
Supergiants are marked as circles, giants as triangles. 
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Fig. 7. — Ratio Voo/vesc as a function of Teg from the current SMC study [solid circles) 
compared with those from the Galactic sample of Lamers et al. (1995, open circles). The 
scalings of Kudritzki & Puis (2000) are also shown [dotted lines). 
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Fig. 8. — Ratio Voo/vgsc as a function of Teff from the current SMC study (solid circles) 
compared with those from the Galactic samples of Kudritzki et al. (1999, open triangles), 
Herrero et al. (2002, open squares) and Repolust et al. (2004, open circles). The scahngs of 
Kudritzki & Puis (2000) are again shown (dotted lines). 



